Magnetic Field Induced Charging Effects in Josephson Junction Arrays 
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A magnetic field induced electric polarization and the corresponding change of an effective junction capacitance 
are considered within a 3D model of disordered Josephson junction arrays. At some threshold field (near the 
Josephson network critical field), the effective junction charge and the related capacitance are shown to reach a 
maximum and to change a sign, respectively. A possibility to observe the predicted effects in artificially prepared 
arrays of superconducting grains is discussed. 
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Recently, a variety of different field-induced phenom- 
ena in high-T c superconductors (HTS), related to their 
extrinsic or intrinsic granularity, have been reported (see, 
e.g., |l|-^| and further references therein). In particular, 
an appearance of an electric-field induced non-zero mag- 
netization signal in granular superconductors (due to the 
Dzyaloshinski-Moria type interaction between an applied 
electric field and an effective magnetic field of circulat- 
ing Josephson currents) was predicted jfj). At the same 
time, an artificially prepared islands of superconducting 
grains, well-described by the various models of Josephson 
junction arrays (JJA), proved useful in studying (both 
theoretically and experimentally) the charging effects in 
these systems, ranging from Coulomb blockade of Cooper 
pair tunneling and Bloch oscillations to propagation 
of quantum ballistic vortices || . 

The present paper addresses yet another related phe- 
nomenon which is actually dual to the above-mentioned 
analog of magnetoelectric effect (described in Ref. ||). 
Specifically, we discuss a possible appearance of a non- 
zero electric polarization and the related change of the 
charge balance in the system of weakly-coupled super- 
conducting junctions (modelled by the random 313 J JAs) 
under the influence of an applied magnetic field. 

To achieve our goal, we apply the so-called model of 
disordered 3D JJAs based on the well-known tunneling 
Hamiltonian (see, e.g., fiof) 
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which describes an interaction between N superconduct- 
ing grains (with phases <pi), arranged in a random three- 
dimensional (3D) lattice with coordinates f\ = (xi, iji, zi). 
The grains are separated by insulating boundaries pro- 
ducing Josephson coupling with energy = J. The 
system is under the influence of a frustrating applied 
magnetic field H . 

Since the field-induced effects considered in the present 
paper are expected to manifest themselves in high enough 
applied magnetic fields (with a nearly homogeneous dis- 
tribution of magnetic flux along the junctions) and as 
long as the Josephson penetration length Xj exceeds the 
characteristic size of the Josephson network d (which is 
related to the projected junction area S, where the field 
penetration actually occurs, as follows S = 7rcP), the 
Josephson current-induced "self-field" effects (which are 
important for a large-size junctions and/or small applied 
magnetic fields) may be safely neglected |T^]. Besides, 
it is known jf2| that in discrete JJAs pinning (by a sin- 
gle junction) actually concurs with the "self-field" effects. 
Specifically, it was found Q that the ratio d/Xj is re- 
lated to the dimensionless pinning strength parameter j3 
as d/Xj — \f(3 suggesting that a weak pinning regime 
(with /J « 1) simultaneously implies a smallness of the 
"self- field" effect and vice versa. And since artificially 
prepared Josephson networks allow for more flexibility 
in varying the experimentally-controlled parameters, it 
is always possible to keep the both above effects down by 
appropriately tuning the ratio d/Xj. Typically [Q-^), in 
this kind of experiments S = 0.01 — 0.1/zm 2 and d <C Xj. 

In what follows, we are interested in the magnetic field 
induced behavior of the electric polarization in a 313 JJA 
at zero temperature. Recall that a conventional (zero- 
field) pair polarization operator within the model under 
discussion reads [Ol 



N 

E 



(6) 



where q L = — 2en^ with m the pair number operator, and 
Ti is the coordinate of the center of the grain. 

In view of Eqs.(l)-(6), and taking into account a usual 
"phase-number" commutation relation, [</>j,rij] = iSij, 
the evolution of the pair polarization operator is deter- 
mined via the equation of motion 



f 



dp 1 r-.„,i 2e A , , 



(7) 



Resolving the above equation, we arrive at the following 
net value of the magnetic-field induced polarization (per 
grain) 



P(H) 



1 

N 

2eJ 
TitN 



< P(t) > 



dt 



' N 



< smfaj^Hyij >, (8) 



where < ... > denotes a configurational averaging over 
the grain positions, while the bar means a temporal av- 
eraging (with a characteristic time r). To consider a 
field-induced polarization only, we assume that in a zero 
magnetic field, P = 0, implying 4>ij(0) = 0. 

To obtain an explicit form of the field dependence of 
polarization, let us consider a site-type positional dis- 
order allowing for weak displacements of the grain sites 
from their positions of the original 31? lattice, i.e., within 
a radius d (which is of the order of the Josephson lat- 
tice parameter and the characteristic junction size, see 
above) the new position is chosen randomly according to 
the normalized distribution function f(f,R). It can be 
shown that a particular choice of the distribution func- 
tion will not change the main qualitative results of the 
present paper. So, assuming, for simplicity, the following 
distribution law 



f(f,R)=f r (f)f R (R), 



with 



and 



(27rd 2 ) 3 / 2 



-(x 2 +y 2 +z 2 )/2d 2 



fn(R) — S(X — d)S(Y — d)6(Z — d), 



(9) 
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we observe that the magnetic field H — (0, 0, H z ) (ap- 
plied along the z-axis) will induce a non- vanishing longi- 
tudinal (along x-axis) electric polarization 



with 



P X (H Z ) = P G(H Z /H ), 



G{z) 



(12) 



(13) 



Here P = edrJ/h, Hq = $o/S with S = ird 2 being 
an average projected area of a single junction, and z — 
H z /H . 

For small applied fields (z <C 1), the induced polar- 
ization P X (H Z ) fts aoH z , where ao = erdJ /HHq is the 
so-called linear magnetoelectric coefficient [Q. However, 
as we mentioned in the very beginning, to correctly de- 
scribe any induced effects in very small external fields, 



both Josephson junction pinning and "self- field" effects 
have to be taken into account ftL2\. 
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FIG. 1. The behavior of the induced effective charge Q/Qo 
(a) and the flux capacitance C/Co (b) in applied magnetic 
field H z /H . 



2 



At the same time, in view of Eq.(6) the induced po- 
larization is related to the corresponding change of the 
effective charge Q = (1/d) J2i < 1i x i > m applied mag- 
netic field as follows 



Q(H Z ) = 



P X (H Z ) 



= Q G(H Z /H ), 



(14) 
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where Qq = erJ/h. 

It is of interest also to consider the related field behav- 
ior of the effective flux capacitance C = TdQ(H z )/d& 
which in view of Eq.(14) reads 



o//j = r (1 (i-2|f 



(15) 



where $ = SH Z , and C = tQ /§ . Figure 1 shows 
the behavior of the induced effective charge Q/Qo (a) 
and the corresponding capacitance C/Cg (b) in applied 
magnetic field H z /Hq. As is seen, at H z /Hq ss 0.8 the 
effective charge reaches its maximum while the capaci- 
tance changes its sign at this field, suggesting a signif- 
icant redistribution of the junction charge balance in a 
model system under the influence of an applied magnetic 
field, near the Josephson critical field Hq. Note that a 
somewhat similar behavior of the magnetic field induced 
charge (and related capacitance) has been observed in 2D 
electron systems fl5j| . 

Taking r = lO^^s for the Josephson relaxation time 
(which is related to the Josephson plasma frequency 
ujp ~ r , known to be the characteristic frequency of the 
system for E c <C J regime, with E c being the Coulomb 
grain's charge energy; typically lu p — 10 9 — 10 11 Hz) 
and J/ks = 90AT for a zero-temperature Josephson en- 
ergy in YBCO ceramics, we arrive at the following esti- 
mates of the effective charge Qo « 10 _16 C, flux capaci- 
tance Co ~ 10~ 11 F, the equivalent current Iq ~ Qq/t 
1(T 6 A, and voltage V ~ Qo/Co ~ 10~ 5 V. We note that 
the above set of estimates fall into the range of param- 
eters used in typical experiments to study the charging 
effects both in single JJs (with the working frequency 
from RF range of u) ~ 10GHz used to stimulate the sys- 
tem H ) and J JAs |l(| suggesting thus quite an optimistic 
possibility to observe the above-predicted field induced 
effects experimentally, using a specially prepared system 
of arrays of superconducting grains. 

In summary a zero-temperature behavior of the in- 
duced electric polarization, effective charge and concomi- 
tant flux capacitance in an applied magnetic field have 
been considered within a 31? model of Josephson junc- 
tion arrays. A possibility of the experimental observation 
of the discussed effects in artificially prepared system of 
superconducting junctions was suggested. 
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